Chemical and Strontium Isotopic Compositions of the Hanjiang Basin Rivers in China: Anthropogenic Impacts and Chemical Weathering   总被引：1，自引：0，他引：1  

Zhifang Xu  Chao Shi  Yang Tang  Hongyin Han 《Aquatic Geochemistry》2011,17(3):243-264

The Hanjiang River, the largest tributaries of the Changjiang (Yangtze) River, is the water source area of the Middle Route of China’s South-to-North Water Transfer Project. The chemical and strontium isotopic compositions of the river waters are determined with the main purpose of understanding the contribution of chemical weathering processes and anthropogenic inputs on river solutes, as well as the associated CO₂ consumption in the carbonate-dominated basin. The major ion compositions of the Hanjiang River waters are characterized by the dominance of Ca²⁺ and HCO₃ ⁻, followed by Mg²⁺ and SO₄ ²⁻. The increase in TDS and major anions (Cl⁻, NO₃ ⁻, and SO₄ ²⁻) concentrations from upstream to downstream is ascribed to both extensive influences from agriculture and domestic activities over the Hanjiang basin. The chemical and Sr isotopic analyses indicate that three major weathering sources (dolomite, limestone, and silicates) contribute to the total dissolved loads. The contributions of the different end-members to the dissolved load are calculated with the mass balance approach. The calculated results show that the dissolved load is dominated by carbonates weathering, the contribution of which accounts for about 79.4% for the Hanjiang River. The silicate weathering and anthropogenic contributions are approximately 12.3 and 6.87%, respectively. The total TDS fluxes from chemical weathering calculated for the water source area (the upper Hanjiang basin) and the whole Hanjiang basin are approximately 3.8 × 10⁶ and 6.1 × 10⁶ ton/year, respectively. The total chemical weathering (carbonate and silicate) rate for the Hanjiang basin is approximately 38.5 ton/km²/year or 18.6 mm/k year, which is higher than global mean values. The fluxes of CO₂ consumption by carbonate and silicate weathering are estimated to be 56.4 × 10⁹ and 12.9 × 10⁹ mol/year, respectively.  相似文献   

Observation targeting with a second-order adjoint method for increased predictability     

Humberto C. Godinez  Dacian N. Daescu 《Computational Geosciences》2011,15(3):477-488

The efficiency of current adjoint-based observations targeting strategies in variational data assimilation is closely determined by the underlying assumption of a linear propagation of initial condition errors into the model forecasts. A novel targeting strategy is proposed in the context of four-dimensional variational data assimilation (4D-Var) to account for nonlinear error growth as the forecast lead time increases. A quadratic error growth model is shown to maintain the accuracy in tracking the nonlinear evolution of initial condition perturbations, as compared to the first-order approximation. A second-order adjoint model is used to provide the derivative information that is necessary in the higher-order Taylor series approximation. The observation targeting approach relies on the dominant eigenvectors of the Hessian matrix associated with a specific forecast error aspect as an indicator of the directions of largest quadratic error growth. A comparative qualitative analysis between observation targeting based on first- and second-order adjoint information is presented in idealized 4D-Var experiments with a two-dimensional global shallow-water model. The results indicate that accounting for the quadratic error growth in the targeting strategy is of particular benefit as the forecast lead time increases.  相似文献   

Analytical model for screening potential CO<Subscript>2</Subscript> repositories     

Roland T. Okwen  Mark T. Stewart  Jeffrey A. Cunningham 《Computational Geosciences》2011,15(4):755-770

Assessing potential repositories for geologic sequestration of carbon dioxide using numerical models can be complicated, costly, and time-consuming, especially when faced with the challenge of selecting a repository from a multitude of potential repositories. This paper presents a set of simple analytical equations (model), based on the work of previous researchers, that could be used to evaluate the suitability of candidate repositories for subsurface sequestration of carbon dioxide. We considered the injection of carbon dioxide at a constant rate into a confined saline aquifer via a fully perforated vertical injection well. The validity of the analytical model was assessed via comparison with the TOUGH2 numerical model. The metrics used in comparing the two models include (1) spatial variations in formation pressure and (2) vertically integrated brine saturation profile. The analytical model and TOUGH2 show excellent agreement in their results when similar input conditions and assumptions are applied in both. The analytical model neglects capillary pressure and the pressure dependence of fluid properties. However, simulations in TOUGH2 indicate that little error is introduced by these simplifications. Sensitivity studies indicate that the agreement between the analytical model and TOUGH2 depends strongly on (1) the residual brine saturation, (2) the difference in density between carbon dioxide and resident brine (buoyancy), and (3) the relationship between relative permeability and brine saturation. The results achieved suggest that the analytical model is valid when the relationship between relative permeability and brine saturation is linear or quasi-linear and when the irreducible saturation of brine is zero or very small.  相似文献   

Updating Markov chain models using the ensemble Kalman filter   总被引：1，自引：0，他引：1  

Dean S. Oliver  Yan Chen  Geir Nævdal 《Computational Geosciences》2011,15(2):325-344

  相似文献   

Three-phase numerical model of water migration in partially frozen geological media: model formulation,validation, and applications   总被引：3，自引：0，他引：3  

Scott L. Painter 《Computational Geosciences》2011,15(1):69-85

Water in the subsurface of the Earth’s cold regions—and possibly the subsurface of Mars—resides in the liquid, vapor, and ice phases. However, relatively few simulations addressing full three-phase, nonisothermal water dynamics in below-freezing porous media have been undertaken. This paper presents a nonisothermal, three-phase approach to modeling water migration in partially frozen porous media. Conservation equations for water (as ice, liquid, and vapor) and a single gas species (in the gas phase and dissolved in water) are coupled to a heat transport equation and solved by a finite-volume method with fully implicit time stepping. Particular attention is given to the method of spatial differencing when the pore space is partially filled with ice. The numerical model is able to reproduce freezing-induced water redistribution observed in laboratory experiments. Simulations of Earth permafrost dynamics and of the formation and evolution of a planetary-scale cryosphere on Mars demonstrate the new capabilities.  相似文献   

G23FM: a tool for meshing complex geological media     

Hussein Mustapha 《Computational Geosciences》2011,15(3):385-397

We present G23FM, a mesh generation tool for discretizing two- and three-dimensional complex fractured geological media. G23FM includes different techniques to generate finite element grids that maintain the geometric integrity of input surfaces, and geologic data and produce optimal triangular/tetrahedral grids for flow and transport simulations. G23FM generates grid for two-dimensional cross-sections, represents faults and fractures, for three-dimensional fractured media, and has the capability of including finer grids. Different examples are presented to illustrate some of the main features of G23FM.  相似文献   

Modelling thermal transients from magmatic underplating—an example from the Vøring margin (NE Atlantic)     

Magnus Wangen  Rolf Mjelde  Jan Inge Faleide 《Computational Geosciences》2011,15(4):771-788

The thermal impact of several kilometre-thick magmatic underplating in the lower continental crust is studied with analytical and numerical methods. Simple analytical solutions are derived for the thermal transient in the case of an infinite depth below the underplate and also for the case of a finite depth (down to the asthenosphere). It is shown that these solutions lead to simple approximations for when the transient surface heat flow is at its maximum, what the maximum is, and for how long the transient lasts. Even though these solutions assume that the underplate is emplaced instantaneously, they are useful in the interpretation of underplating over finite time spans. A numerical scheme is suggested for the modelling of underplating that handles both short time intervals as well as long intervals. The scheme treats magmatic underplating in a mass and energy conservative manner, and it is compared against the analytical solutions. Finally, the analytical and numerical results for thermal transients are applied to a transect from the Vøring margin (NE Atlantic), with respect to various degrees of early Cenozoic magmatic intrusion. It appears that more than half of the lower crustal body (LCB) in the Vøring margin must be magmatic underplating for the vitrinite reflectance to be substantially higher than for the non-magmatic case, where the LCB is assumed to comprise Caledonian crust.  相似文献   

Magmato-sedimentary Carboniferous to Jurassic evolution of the western Tauern window,Eastern Alps (constraints from U-Pb zircon dating and geochemistry)     

Veselá Petra  Söllner Frank  Finger Friedrich  Gerdes Axel 《International Journal of Earth Sciences》2011,100(5):993-1027

New geochronological U-Pb (LA-ICP-MS) zircon data and geochemical analyses from the Variscan orthogneisses and metavolcanic rocks in the western Tauern window are presented and used to reconstruct the pre-Alpine evolution of this area. The late- and post-Variscan stage in the Tauern window was characterised by distinct magmatic pulses accompanied by the formation of volcano-sedimentary basins. The magmatic activity started in the Visean (335.4 ± 1.5 Ma) with the intrusion of a K-rich, durbachitic biotite-granite (protolith of the Ahorn gneiss). Following a period of exhumation and erosion, Westfalian–Stefanian volcanics were deposited (Grierkar meta-rhyodacite: 309.8 ± 1.5 Ma; Venntal meta-rhyolite: 304.0 ± 3.0 Ma). A renewed magmatic pulse occurred in the Early Permian, producing large volumes of tonalites and granodiorites (Tux meta-granodiorite: 292.1 ± 1.9 Ma). The youngest magmatism is characterised by pyroclastic and tuffitic deposits (Pfitsch meta-rhyolite: 280.5 ± 2.6 Ma; Schönach valley meta-andesite: 279.0 ± 4.8 Ma). This volcanism was probably related to crustal extensional faulting within an intra-continental graben and horst setting, asthenospheric upwelling and heat flow increase due to the onset of the Permian rifting. The Permo-Triassic peneplanation and subsidence is documented by shallow marine and evaporitic deposits. Probably in the Middle Jurassic times, the area was flooded and in the Late Jurassic the whole area was covered by limestones, representing post-rift sediments on the southern European continental margin.  相似文献   

Geochemical evaluation of Niger Delta sedimentary organic rocks: a new insight     

Akinsehinwa Akinlua  Nelson Torto 《International Journal of Earth Sciences》2011,100(6):1401-1411

A geochemical evaluation of Niger Delta organic matter was carried out using supercritical fluid extraction (SFE) sample preparation procedure. Comparison of geochemical significance of gas chromatographic data of rock extracts of SFE with those of Soxhlet extraction method from previous studies was made in order to establish the usefulness of SFE in geochemical exploration. The assessment of geochemical character of the rock samples from the comparison and interpretation of other geochemical parameters were used to give more insights into understanding the source rocks characteristics of onshore and shelf portions of the Niger Delta Basin. The results of the gas chromatographic (GC) analysis of the rock extracts across the lithostratigraphic units show that Pr/Ph, Pr/nC₁₇, Pr/nC₁₈, CPI and odd/even preference ranged from 0.07 to 12.39, 0.04 to 6.66, 0.05 to 13.80, 0.12 to 8.4 and 0.06 to 8.12, respectively. The Rock-Eval pyrolysis data and geochemical ratios and parameters calculated from the GC data showed that most of the samples are mature and have strong terrestrial provenance while a few samples have strong marine provenance. The few marine source rocks are located in the deeper depth horizon. Pr/Ph and standard geochemical plots indicate that most of samples were derived from organic matter deposited in less reducing conditions, i.e. more of oxidizing conditions while a few samples have predominantly influence of reducing conditions. The results of trace metal analysis of older samples from Agbada Formation also indicate marine and mixed organic matter input deposited in less reducing conditions. The results obtained in this study are comparable with those obtained from previous studies when Soxhlet extraction method was used and also indicated the presence of more than one petroleum systems in the Niger Delta.  相似文献   

Lattice-preferred orientations of late-Variscan granitoids derived from neutron diffraction data: implications for magma emplacement mechanisms     

Axel Müller  Bernd Leiss  Klaus Ullemeyer  Karel Breiter 《International Journal of Earth Sciences》2011,100(7):1515-1532

The lattice-preferred orientation (LPOs) of two late-Variscan granitoids, the Meissen monzonite and the Podlesí dyke granite, were determined from high-resolution time-of-flight neutron diffraction patterns gained at the diffractometer SKAT in Dubna, Russia. The results demonstrate that the method is suitable for the LPO analysis of polyphase, relatively coarse-grained (0.1–6 mm) rocks. The Meissen monzonite has a prominent shape-preferred orientation (SPO) of the non-equidimensional minerals feldspar, mica and amphibole, whereas SPO of the Podlesí granite is unapparent at the hand-specimen scale. The neutron diffraction data revealed distinct LPOs in both granitoids. The LPO of the non-equidimensional minerals feldspar, mica and amphibole developed mainly during magmatic flow. In the case of the Meissen monzonite, the magmatic flow was superimposed by regional shear tectonics, which, however, had no significant effect on the LPOs. In both samples, quartz shows a weak but distinct LPO, which is atypical for plastic deformation and different in the syn-kinematic Meissen monzonite and the post-kinematic Podlesí granite. We suggest that, first of all, the quartz LPO of the Meissen monzonite is the result of oriented growth in an anisotropic stress field. The quartz LPO of the Podlesí granite, which more or less resembles a deformational LPO in the flattening field of the local strain field, developed during magmatic flow, whereby the rhombohedral faces of the quartz crystals adhered to the (010) faces of aligned albite and to the (001) faces of zinnwaldite. Due to shape anisotropy of their attachments, the quartz crystals were passively aligned by magmatic flow. Thus, magmatic flow and oriented crystal growth are the major LPO-forming processes in both granitoids. For the Meissen monzonite, the solid-state flow was too weak to cause significant crystallographic re-orientation of the minerals aligned by magmatic flow. Finally, the significance of our results for the evaluation of the regional tectonic environment during magma emplacement is discussed. The discussion on the regional implications of the more methodologically oriented results provides the basis for future, more regionally aimed studies in view of the fabric characteristics of such plutons and their developing mechanisms.  相似文献